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Equilibrium director structures in two thin hybrid planar films of biaxial nematics 
are investigated through Markov chain Monte Carlo simulations based on a lattice 
Hamiltonian model within the London dispersion approximation. While the sub¬ 
strates of the two films induce similar anchoring influences on the long axes of the 
liquid crystal molecules (viz. planar orientation at one end and perpendicular, or 
homeotropic, orientations at the other), they differ in their coupling with the minor 
axes of the molecules. In Type - A film the substrates do not interact with the minor 
axes at all (which is experimentally relatively more amenable), while in Type - B, 
the orientations of the molecular axes at the surface layer are influenced as well by 
their biaxial coupling with the surface. Both films exhibit expected bending of the 
director associated with ordering of the molecular long axes due to surface anchoring. 
Simulation results indicate that the Type - A film hosts stable and noise free director 
structures in the biaxial nematic phase of the LC medium, resulting from dominant 
ordering of one of the minor axes in the plane of the substrates. High degree of 
this stable order thus developed could be of practical interest for in-plane switching 
applications with an external field. Type - B film, on the other hand, experiences 
competing interactions among the minor axes, due to incompatible anchoring influ¬ 
ences at the bounding substrates, apparently leading to frustration, and hence to 
noisy equilibrium director structures. 
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INTRODUCTION 


The biaxial nematic phase (Nb) of liquid crystals (LC), predicted theoretical 
very early JlL and realised experimentally in the past decade, in bent-core w. 
tetrapode ^7|, and polymeric js, 9| systems, is characterised by a primary director 
n and a secondary director m (perpendicular to n). Field-induced switching of the 
secondary director is envisaged to be faster than the primary director in the biaxial 
nematic phase, a fact which endows these nematics with a promising potential for 
use in fast switching electro-optic devices 10|, ]T|. The orthorhombic Nb phases with 


D 2 h symmetry are suggested to be desirable 


12| for ready applications. 


While the current experimental studies are still concerned with unambiguous con- 
hrmation of macroscopic biaxiality 3^14|, computer simulations 151 on the other 
hand have been playing a signihcant role in the investigation of biaxial systems. A 
molecular dynamics (MD) simulation of the bulk biaxial Gay-Berne fluid under the 
action of an electric held has convincingly shown that the switching of the sec¬ 
ondary director is an order of magnitude faster than the switching of the primary 
director. 

In the biaxial nematic phase, a different pathway for fast switching between dif¬ 
ferent birefringent states (compared to conventional uniaxial LC systems) is possible 
because the birefringence can be changed by a rotation of the short axes which are 
thermally ordered, while the orientation of the long axes could be kept hxed [l^. 
A possible device conhguration to achieve this objective is to use a him of biaxial 
liquid crystal conhned in a planar cell with hybrid boundary conditions, wherein the 
geometry could constrain the orientation of the primary director (of the long axes) in 
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the biaxial nematic phase, leaving the secondary director (of one of the short axes) 
free for switching with an appropriate (in-plane) held. Studies on uniaxial hybrid 


hlms have established 


18 


19| that a bent-director conhguration could be realised if 


the him thickness is greater than a critical thickness determined by the curvature 
elasticity of the medium and the surface interaction strength. Preliminary work on 
their biaxial counterpart was also carried out earlier 20|. 

In this context, we investigated the equilibrium director structures in two planar 
hlms of biaxial liquid crystals, in the uniaxial and biaxial phases of the medium. 
The anchoring inhuences of the two substrates comprising the cell are used to pin 
the orientation of the primary director (ordering direction of the major (long) axes 
of the molecules) near the two surfaces so as to result in a bent-director hybrid 
structure. The substrates can be chosen either not to have inhuence on the minor 
axes of the molecules (pure uniaxial coupling with the substrate, say Type - A him), 
or to couple with the minor axes as well (biaxial coupling with the substrate, say 
Type - B him). We simulated director structures in both these hlms based on Markov 
chain Monte Carlo (MC) sampling technique, constructing corresponding equilibrium 
ensembles. This paper reports our results, examining the role of biaxial anchoring on 
the orientational ordering in the medium, and looking for their suitability to utilize 
the secondary director for possible switching applications. 

In section II we introduce the lattice model of the medium, and details of the 
anchoring conditions of the two hlms. The MC simulation are also discussed in this 
section. The equilibrium director structures of the him in the two nematic phases 
obtained from the computations are depicted and discussed in section III. We also 
examine the ehect of varying the cell thickness, as well as of the relative anchoring 
strengths at the two substrates on the director structures in these hlms. The last 
section summarizes our conclusions. 
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II. MODEL AND SIMULATION DETAILS 


We consider a planar hybrid film comprising of LC molecules with D 2 h symme¬ 
try. We assign the right-handed triad {X, y,Z} to represent the laboratory-fixed 
frame and {x, y, z} to represent the molecular-fixed frame. We let 2 : direction rep¬ 
resent the molecular long axis, while the other two (minor) axes are represented by 
y and x. The film is obtained by confining the biaxial molecules between two planar 
substrates taken to be in the X-Y plane. Fig. [1] shows the schematics of a biaxial 
molecule, orienting influences at the two substrates of the planar film, and the ref¬ 
erence axes of the laboratory. The orientational interactions between LC molecules, 
relevant to the present study, are conveniently accounted for, by adopting a lattice 
Hamiltonian model wherein the molecules are located at the lattice sites represented 
by unit vectors {x, y, z} specifying the individual molecular orientations. Within this 
lattice description, in a film of thickness d the substrate planes are positioned at 
Z = 0 and Z = d + 1 (lattice units). The influence of the substrates is simulated 
by introducing two bounding layers of molecules contained in these planes with the 
designated, but fixed, orientations, referred to in the literature as ghost molecules 


2l[ |. The anchoring conditions in both the films are such that the long axes of the 


molecules are hybrid-aligned (planar orientation parallel to say, y-axis at one sub¬ 
strate and homeotropic at the other, parallel to z-axis). For adequately thick films, 
the primary director n is bent satisfying the two incompatible boundary conditions. 
We distinguish two scenarios: (a) in Type - A film, the ghost molecules interact with 
the LC molecules in the surface layer anchoring only their long axes, thus implying 
that the substrate hosts only cylindrical symmetric rod-like LC constituents, and (b) 
in Type - B film the ghost molecules themselves have D 2 h symmetry and interact 
through a Hamiltonian model appropriate to the biaxial system.. These correspond 
to two qualitatively distinct chemical treatments of the anchoring substrates. 
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FIG. 1: Schematic illustrations of (a) Typical biaxial molecule (b) anchoring directions at 
the two substrates of the film: homeotropic at the lower (k=l) substrate and planar at the 
top (k=d) (c) orientation of the laboratory frame relative to the film. 


In Type - A film, the molecular z-axes are anchored along the Z-direction in the 
surface layer near, say, lower substrate of the cell (k=l), while they are kept planar 
(say, parallel to Y-axis) in other surface layer, near the substrate (k=d). As the 
substrates do not interact with the minor axes of the molecules in this him, they are 
not, a priori, oriented in any specihc direction, until guided by equilibrium conditions 
under a given interaction model. 

In the Type - B him, we impose anchoring conditions of equal strength on all the 
three molecular axes at each of the two substrates, and use the biaxial Hamiltonian 
model to account for their interactions with the substrates. The boundary conditions 
on the him at the two surfaces is summarized as (see Fig. [1]): 
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Substrate 1 at Z=0 ■. z \\ Z,y \\Y,x \\ X 
Substrate 2 at Z=d+1 ■. z \\Y,y \\ Z,x \\ X 

A. Model Hamiltonian 


The biaxial LC molecules are assumed to interact through a pair-wise additive 
lattice Hamiltonian 22|] within the London dispersion approximation, expressed in 
terms of generalised Wigner rotation matrices as: 


Uiuij) — —ejj{P2(cos(/3jj)) + 2 [-Rg2(wp) + -R2o(‘^p)] 

+4 


( 1 ) 


where Cjj = e sets the energy scale, and is used to dehne the reduced temperature, 
u{a,j3,'j) are the set of Euler angles which specify the rotations to be performed 
in order to bring the reference frame of two molecules i and j in coincidence, R^n 
are symmetrized Wigner functions, P2{cos{/3ij)) is the second Legendre polynomial 
and Xd quantihes the biaxial interaction between the molecules. The average values 
of dehne the order parameters of the medium in the nematic phases. These 
are: the uniaxial order < R^q > (along the primary director), the phase biaxiality 
< i? 2 o >, and the molecular contribution to the biaxiality of the medium < R^ o > , 


and the contribution to uniaxial order from the molecular minor axes < Rq 2 > 


22 |. 


For simulation purposes, the above Hamiltonian is recast in the cartesian form , as 


U = -e{-H33 


\/6 Arf (Hn - H22) + A^ (Hn + H22 - W2 - f^ 2 i) - -}• ( 2 ) 
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Here, Vab = {ua-VbY, and the unit vectors Ua,Vb, [a,b = 1,2,3], are the three axes 

of the two interacting neighbouring molecules. sets the relative importance of 

the biaxial interaction in the Hamiltonian, while e (set to unity in the simulations) 

/ ksT* 

dehnes the temperature scale (T =-) , where T* is the laboratory temperature 

in Kelvin. 

This model in a bulk system was studied extensively both through mean-held 
analysis and MC simulations based on Boltzmann sampling methods. The value of 
Xd for the present study is kept at 0.35, keeping in view the high degree of biaxiality 
it induces, as well as the convenience of a wider biaxial nematic range of temperature 
made available for our study [221]. 


B. Simulation Details 

A planar hybrid him of (lattice) dimensions 15 x 15 x d (d = 6, 8, 10, 12 layers) is 
considered in the present work. Periodic boundary conditions are applied along the 
laboratory X and Y directions, so as to minimize hnite size ehects. The anchoring 
conditions applied at the two substrates (contained in the X-Y plane) depend on 
the specihc choice of the him (Type A or B), and their relative strengths are chosen 
as desired. We index the LC layers with k, starting from the substrate imposing 
homeotropic anchoring inhuence on the long molecular axes. The LC ghost molecules 
in the substrate layers (which are adjacent to the two bounding layers of the LC 
medium) do not participate in the Monte Carlo dynamics. The interaction strengths 
of the long molecular axes at the two substrates are represented by ei and e^, and 
ed can be varied relative to ei (which is set equal to l).The temperature T' (in 
dimensionless units) is set by the coupling strength e in the Hamiltonian (EqnJT]). 

In each simulation, the reduced temperature is varied from 2 to 0.05 in steps of 
0.005, and at each temperature the him is equilibrated for 5 x 10^ lattice sweeps 
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FIG. 2: Variation with temperature of system properties a bulk biaxial film (15 x 15 x 15): 
specific heat Cy peaks mark the phase transitions I — Njj and Njj — Nb- Uniform growth 
of the order parameters major order parameters i?Qg, R 22 can be observed. 


(attempted moves over all the sites), while the data from the resulting canonical 
ensembles are collected over a production run of 5 x 10® lattice sweeps. The physical 
properties computed are the average values of energy E, the specihc heat Cv, the 
order parameters RIq, RI 2 , R^o^ ^22 as well as their susceptibilities, and layer- 


wise orientation of the local directors (orientation averaged over the layer) with 
respect to the laboratory Z-axis (polar angle 6, layerwise) as well as the angle made 
by the local layer-wise director with respect to laboratory X-axis (azimuthal angle 
(j), layerwise). We computed both layer-wise properties (to examine the director 
structure and its relative changes) as well as the him properties (averaged over the 
sample), as a function of temperature, for a hxed relative anchoring strength and 
layer thickness. The simulations are then repeated by varying the thickness effecting 
the length scale of the system, and also varying the anchoring strength (e^) to look 
for possible anchoring-induced transitions among the different director structures. 
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III. RESULTS AND DISCUSSION 


Fig. [2] depicts the variation of specific heat and the major order parameters in a 
bulk biaxial fluid (without conhnement) of dimensions (15 x 15 x 15) computed based 
on MC simulations under periodic boundary conditions, using the above dispersion 
Hamiltonian model at Ad=0.35. The two peaks in the specihc heat at Ti = 1.123 
and T 2 = 0.571 indicate the high temperature / — Nu and low temperature Nu — Nb 
phase transitions 22|]. The phases, as identihed by the Metropolis-based Boltzmann 
sampling- procedure employed here, are marked by uniform growth of the uniaxial 
order R^q and the biaxial order R 22 , as a function of temperature, bringing out the 
onset of the uniaxial phase Nu and low temperature biaxial phase Nb- Fig. E] shows 
their variation in the Type - A him. Fig. H] provides variation of all the four order 
parameters and their susceptibilities in Type - A him, averaged over the sample. 
The qualitative diherence between the thermal evolution of these order parameters 
in the bulk and Type - A him (in particular Rqq) may be noted for further discussion 
later. The director structure in the him is better appreciated by examining the 
degree of ordering in each layer and the orientation of the local layer director with 
respect to laboratory axes, specihed by their 9 and 0. These layer-wise properties 
are shown for Film-A in Fig. [5l Similar data were obtained for Film-B as well, and 
are represented in Figs. [6]- IH respectively. It is to be noted that layer-wise angles 
of the local directors are obviously not meaningful in the isotropic phase, and hence 
such information in these hgures (Figs. [5] and E]) (generated automatically during 
computations) are to be ignored. Such data are relevant only in the ordered phases 
and are subject of our discussion. 



10 



FIG. 3: Variation with temperature of system properties in Type - A film. Specific heat 
Cv peaks mark the phase transitions I — Njj and Njj — Nb- The bent - director structure 
of uniaxial order parameter Rqq in the uniaxial phase and the increase in biaxial order R 22 
in the uniaxial phase itself, are attributed to effects of geometric confinement. 

A. Type - A Film 


We note at the outset that in this film the substrate interacts with only the long 
axes of the molecules, while LC molecules themselves have biaxial interaction among 
them. This system undergoes two transitions at Ti = 1.021, and the second at 
T 2 = 0.525, exhibiting the sequence Nb — Nu — I (see Fig. [3]). It is observed that 
both transition temperatures are lowered, in comparison with the those in a bulk 
sample (see Fig. [2l). 


Figs. 4(a) and 4(b) show the four order parameters of the system along with 
their susceptibilities. It is observed from Fig. 4(a) that the susceptibility prohle, 
which is a measure of the fluctuations in the long-range order of the primary director, 
shows two peaks. The larger peak (at Ti) corresponds to rapid increase of i?QQ in 
the uniaxial nematic phase. The R^q curve displays a bent-director structure and 
attains a maximum value of 0.6 in this phase. The smaller peak close to T 2 signals 
a sudden change in the slope of the Rqq curve leading to a steady increase of the 
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FIG. 4 : Variation with temperature of the four system order parameters (lines with sym¬ 
bols) and their susceptibilities (continuous lines) in Type - A film (a) Rqq and R22 (b) Rq2 
and R2q. 


dominant (uniaxial) order towards a maximum value of unity deep in the biaxial 
phase. The susceptibility of the biaxial order, on the other hand, shows a single 
broad peak at a temperature slightly higher than T 2 , signalling the onset of biaxial 
order in this phase. The non-zero value of R 22 at the onset of the I — Nu transition 
itself is attributed as due to the geometric conhnement (cf. Fig. |2]of bulk system). 

Contribution to the ordering along the uniaxial director originating from molecular 
minor axes, i?Q 2 and the phase biaxiality parameter arising from the non-cylindrical 
distribution of the molecular long axes around the primary director, R\q are shown 


in Fig. 4(b) As compared to their values in a bulk system, we observe that these 
parameters have non-zero values in the uniaxial nematic phase, - a clear manifestation 
of the geometric conhnement. 

The layer-wise variations (for the layer index A: = 1 to d) of properties connected 
with director structures, plotted along with the corresponding bulk values of the 


him (sample averages) are shown in Figs. 5(a) - 5(c) It is observed from Fig. 5(a) 
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FIG. 5 : Layer-wise variation of different properties of the Type A film, with temperature: 
(a) i?gQ and i?22 (b) R02 i?20 (^) Angle 9 and Angle cj) (see text). 
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that layer-wise Rqq values vary smoothly, unlike the abrupt jump at T 2 exhibited 
by the sample average. Further, the order values in all the layers asymptotically 
reach the maximum value of 1.0 at the lowest temperature. The variation also shows 
that the middle layers (A: = 4 and 5), being the least influenced by the substrate 
boundaries, are most effective in contributing to the critical onset of the order at the 
transition. The sample average and layer-wise behaviour of biaxial order shows that 
in the middle layers R 22 starts increasing from the I — Njj transition itself, but a 
signihcant increase is observed only at the Nu — Nb transition. The layer-wise R ^2 


and i? 2 o are shown in Fig. 5(b) It is observed that they have small, but non-zero 
values in the uniaxial nematic phase in the middle layers of the him. 

The plots of the layer-wise 6 (angle made by the layer-wise director with respect 
to the laboratory Z-axis) and layer-wise cp (angle made by the layer-wise director 


with respect to the laboratory X-axis) are shown in Fig. 5(c), Focusing on the data 
only below the clearing point, it is observed that as the temperature is lowered from 
the isotropic phase, layer-wise 9 values stabilise at certain hxed values in the nematic 
phase. The value of this angle increases monotonically from k = 1 layer to the k = d 
layer, clearly indicating that the layer-wise primary director bends gradually from the 
lower substrate to the upper substrate, corresponding to the expected bent-director 
structure. The layer-wise (p angles show a sudden flip by approximately 90° just 
below the / — Nu transition, particularly in the region of middle layers, pointing 
to the onset of a long-range order, essentially conhning the primary director to the 
laboratory YZ plane, as the system is cooled. The prohle of the layer-wise angle also 
suggests that this bent-director structure is not affected by the Nu — Nb transition. 
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FIG. 6 : Variation with temperature of system properties in Type-B film - peaks of the 
specific heat Cy mark the phase transitions I — Njj and Nu — Nb- The phases are identified 
by the growth of the major order parameters i?Qg and R22- 

B. Type-B Film 


Introduction of biaxial coupling between the substrates and the surface him layers, 
in the presence of an already constrained bent primary director, introduces incom¬ 
patible boundary conditions on the minor axes as well, and seem to lead to some 
curious manifestations. 

Referring to Fig. |6l it may be noted that the onset of biaxial order in Type-B him 
occurs only at the Nb — Nb transition unlike the case of Type - A him. Further, 


Figs. 7(a) and 7(b) show that all order parameter susceptibilities, except that of R 22 , 
display appreciable huctuations in the biaxial phase, and these are rehected on their 
(average) order prohles as well. 

The system and layer-wise variations of diherent order parameters are shown in 


Figs. 8(a) - 8(b) The layer-wise Rqo values in Fig. 8(a) do not huctuate in the 


biaxial phase and attain maximum ordering, interestingly unlike the system order 
parameter. The layer-wise biaxial order shows a non - zero value just below the 
/ — Nu transition and grows gradually, whereas the system biaxial order develops 



















15 



FIG. 7 : Variation with temperature of the four system order parameters (lines with sym¬ 
bols) and their susceptibilities (continuous lines) in Type-B film (a) i?Qg and i?22 (b) -^02 
and i?20- 


appreciably only at the Nu — Nb transition. Both fluctuate significantly in the biaxial 
phase similar to the other order parameters. The layer-wise i ?02 and R\q shown in 


Fig. 8(b) have small, nonzero values in the uniaxial nematic phase, but show similar 
fluctuations in the biaxial phase. 

The layer-wise 9 and layer-wise (p values shown in Figs. |8(cJ] depict variations with 
temperature in the two nematic phases, which are largely similar to the behaviour of 
Type - A him, but for signihcant huctuations on the onset of the Nu — Nb transition. 
It appears that the layer-wise magnitudes of the two dominant orders are relatively 
stable, but their orientations are not. 

We now discuss the evolution of the uniaxial sample order (i?oo) fh® two dims as 
the system transits through the two transitions, specihcally in comparison with the 


layerwise behaviour (Figs. 5(a), 8(a)). Focussing on Type - A him initially, we note 


that the sharp increase in this order at the second Nu — Nb transition is unlike the 
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FIG. 8 : Layer-wise variation of different order parameters of the Type-B film, with tem¬ 
perature: (a) TJqq and i?22 (b) R02 i?20 (^) Angles 9 and 4 > (see text). 
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Mont Carlo steps 


(a) (b) 

FIG. 9: Type - A film (a) Variation of maximum eigen values (qmax) associated with each 
of the ordering tensors of the three molecular axes (x, y, z), as a function of temperature 
(b) The orientations of the corresponding eigen vectors plotted as a function of the Monte 
Carlo steps after equilibration, at temperature T = 0.5 below the Njj — Nb transition 
temperature. 

variation in the bulk sample (Fig. [2]), and is not supported by its layer-wise variations. 
The latter are smoothly varying across the transition, and the director angles {6, 0) 
also do not betray the presence of an intermediate phase transition. And yet, the 
bulk order of the him, dehned and computed as corresponding to the maximum 
eigen value of the ordering tensors of the three molecular axes, is very sensitive to 
this transition. The case of the Type - B him is qualitatively the same, but for the 
onset of signihcant huctuations in the Nb phase at T 2 . This points to the need to 
further investigate this system in terms of thermal variations of the eigen values of 
the three axes separately and examine their behavior across this transition. 

Accordingly, we computed equilibrated averages (over the him) of the maximum 
eigen values Qy, Qz) of the ordering tensors {Qxx, Qyy, Qzz) of the three molecular 
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FIG. 10: Type - B film ;(a)Variation of maximum eigen values (qmax) associated with each 
of the ordering tensors of the three molecular axes (x, y, z), as a function of temperature; 
(b) The orientations of the corresponding eigen vectors plotted as a function of the Monte 
Carlo steps after equilibration, at temperature T = 0.5 below the Njj — Nb transition 
temperature. 


axes and depict their variation for the case of Type - A film in Fig. 9(a) We also 
show the directions {9, (p) of the corresponding three eigen vectors, as a fnnction 


of MC steps after eqnilibration, in the Nb phase, in Fig. 9(b) From these two 
figures, it is evident that the onset of the biaxial phase in this confined system (at 
r2=0.5) leads to a sudden change in the direction of dominant order of the film 
itself. While the alignment of the long molecular axes defines the primary director 
till T 2 , it is the molecular y axes which are the most ordered among the three, below 


this temperature. In conjunction with Fig. 9(b), we see that the ordering direction 


of this calamitic axis is indeed in the laboratory X-direction in the biaxial phase, 
while the other two eigen vectors are confined to the laboratory YZ plane, mutually 
perpendicular to each other. It may be noted that the onset of a biaxial phase thus 
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leads to maximal ordering of the second major axes, wholly contained in the plane 
of the substrate, and the anchoring conditions imposed in this case constrain only 
the long axes of the LC molecules, along with the minor a;-axes of the molecules. 

In contrast. Type - B him which imposes anchoring restrictions on the minor axes 
of molecules as well at both the ends, presents a very different scenario. We refer 


to Fig. 10(a) showing Qy, Qz) as a function of temperature in this him. At the 
onset of the second transition ( at T 2 0.51), Qz ( 0.617) is higher than the Qx 

(~ 0.536) and Qy (~ 0.447). However on further cooling, Qx crosses the value of Qz 
(~ 0.632 at T' = 0.433), while Qy remains less than g^. At lower temperatures Qx 
saturates at ~ 0.87, while Qy and qz saturate to a value of < 0.72. Eigen vector of qz 
makes an angle 6 ~ 45° with the laboratory Z-direction and (p ~ 90° with laboratory 
X-direction, thereby indicating that the bent-director structure originating from the 
ordering of the molecular 2 ;-axes is contained in the YZ plane, as is also the case in 
the high temperature nematic phase. Curiously the eigen vector of qy is oriented at 


angles 9 ~ 45° and with a fluctuating cp varying between 0° to 90° (see Fig. 10(b)). 
The maximal ordering direction in this him is determined by the ordering tensor of 
the molecular x-axes and its azimuthal angle (p is eventually contained in the plane 
of the substrate pointing to the laboratory X-direction. It may be noted that the 
corresponding huctuations of the azimuthal angle (p of the local directors of the x 


and y molecular axes are complementary (Fig. 10(b)). 

It is thus obvious that the observed signihcant noise in the thermal averages of 
i?Q 0 , 7?2o Rq 2 at the onset of the Nb phase arises due to the huctuations in the 


directions of the diherent ordering tensors (Figs. 7(a) and 7(b) The relative stability 
observed in the case of R 22 (in the Nb phase) perhaps follows from the interesting 
circumstance that the biaxial order in this him is dehned by the ordering of the 
molecular y and 2 : axes; while the ordering direction of the 2 : axes is fairly stable, that 
of y axes huctuates however complementary to the primary director dehned by the 
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molecular x axes. The qualitatively different scenario of the Type - B him, relative to 
Type - A, arises due to the imposition of additional restrictions on all the molecular 
axes at the two substrates. 

Comparing the two hlms in their biaxial phase, we observe that both have primary 
director (dehned as the direction of maximum molecular order) contained in the 
plane of the substrates, with molecular y-ax.es dehning such a direction for the Type 
- A him while x-axes play that role for Type - B him. However, imposition of 
anchoring constraints on all molecular axes leaves the second him frustrated, leading 
to signihcant huctuations of the ordering directions. Moreover from a practical point 
of view, it is much easier to prepare substrates which need to restrain only one type 
of axes of the system. Thus Type - A lends itself as a possible structure containing 
stable and signihcant order contained within the plane of the him (and hence can be 
unambiguously coupled to an in-plane external held), with the added convenience in 
the preparation of required substrates. 


C. Effect of thickness in a Type - A film 


Because of its simplicity and stable structures, we examined the ehect of varying 
the thickness 'd in the Type - A him, {d = 6, 8, 10, 12 lattice units), while re¬ 
taining the same lateral dimensions, and anchoring strengths at the two substrates. 


Fig. [TT] shows this ehect on the specihc heat Cy and Figs. 12(a) and 12(b) depict the 
consequent changes in i?QQ, -^02 respectively. 

It is observed that the specihc heat prohles become sharper as the thickness 
increases (size ehect). The uniaxial order parameter i?Qg increases marginally in the 
uniaxial phase as the thickness increases, retaining its bent director structure for all 
thickness. The biaxial order i ?22 shows a slight decrease as the thickness increases 
and shows a marked increase in the temperature of the nematic phase itself for a 
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FIG. 11 : Comparison of variations of the specific heat Cy as a function of temperature 
in films of thickness d (in lattice units) = 6 (stars), 8 (circles), 10 (upward triangles) and 
12 (diamonds), in Type - A film. 
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(a) (b) 

FIG. 12 : Comparison of variation of different order parameters with temperature in Type 
- A film for different film thickness {d). (a)i?QQ and (t*) ^02 d = 6 (stars), d 

= 8 (circles), d = 10 (upward triangles) and d = 12 (diamonds), in Type - A film. 

thinner film {d=6). The Rq 2 and i? 2 o values decrease as the thickness increases (size 
effect). These observations indicate that the gross features of director structures in 
Type - A him are relatively insensitive to the him thickness, alluding to some degree 
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(a) 


(b) 


FIG. 13 : Variation of order parameters with temperature for different anchoring strengths 
(e^z) in Type A film of thickness 6 lattice units: (a) and i?22 (b) Rq2 R\q. 


of flexibility in its design. 


D. Effect of anchoring strength in Type - A film 


We further examined the Type - A him (at hxed thickness d = 6) with respect 
to a change in the anchoring strength at the surface layer. We relax the strong 
anchoring condition (e^ = 1) and vary its value now from 0.0 to 0.6 in steps of 0.1, 
and compute equilibrium averages (as a function of temperature) of the four order 
parameters. We depict their dependence on in Fig. [131 The other simulational 
conditions remain the same as before. 


We observe from Figs. 13(a) and 13(b) that for anchoring strengths 6^=0.0 (no 


anchoring influence at the top substrate) and for a low value of = 0.1, the uniaxial 
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order i?Qg and the biaxial order R 22 attain maximnm valnes of 1.0 and 0.5, respec¬ 
tively, - for example withont the characteristic featnres observed earlier (with = 1) 
at the onset of the Nb phase. Their variations are similar to bnlk LC systems with¬ 
ont conhning snrfaces (see Figs. [2] and [3]). For > 0.2 the primary director assnmes 
a bent strnctnre and the primary order is constrained with an npper bonnd of 0.8 
for Sd = 0.2 in the nniaxal nematic phase. A similar sharp difference is exhibited by 
i ?22 well, and its temperatnre variation qnalitatively changes above the threshold 


valne of = 0.2. The order parameters RI 2 and R^q, shown in Fig. 13(b), start with 
very low valnes at = 0.0 and 0.1, bnt increase signihcantly for higher anchoring 
strength (e^ > 0.2). Thns it appears that a minimnm threshold anchoring strength 
(q > 0.2) is necessary for the him to exhibit the cnrions strnctnres, described ear¬ 
lier, arising from the bent formation of its director. The progressive development of 
this scenario with increase of is evident from the gradnal decrease of the system 


primary order (Fig. 13(a)) in the nniaxial phase. 


IV. CONCLUSIONS 

Eqnilibrinm director strnctnres in thin planar dims of biaxial LC medinm, im¬ 
posing hybrid bonndary conditions at the two snbstrates on the molecnlar long axes 
(z-axes), are investigated throngh Monte Carlo simnlation, based on a Hamilto¬ 
nian model nnder dispersion approximation (at Xd = 0.35). Boltzmann sampling is 
adopted which corresponds to the established phase seqnence involving nniaxial and 
biaxial phases {Nb — Nu — I) in bnlk samples. Two different dims were considered: 
Type - A in which only the long axes are anchored at the snbstrates and Type - B 
in which all molecnlar axes are snitably anchored. Very interestingly. Type - A dim 
leads to a stable, noise free, high order parameter valne for Rqq deep in the biaxial 
phase region, and this order is dne to the near perfect alignment of the molecnlar 
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y-axes in a direction perpendicular to the plane hosting the bent director structure 
associated with the long molecular axes. This scenario is appealing for two reasons: 
the Hamiltonian model is realistic as it accommodates both the perturbing inter¬ 
actions (the last two terms in EqnH]), and secondly the substrates are needed to 
influence suitably only the long molecular axes, which is more practical. Under these 
convenient conditions we hnd that in the biaxial phase, one can generate an in-plane 
order of high degree involving one of the minor axes of the molecules, while the bent 
director from the anchored alignment of the long axes is conveniently locked within 
the plane perpendicular to this dominant order. Type - A him appears to be a good 
candidate for further studies of switching with external in-plane helds. Type - B him 
on the other hand is frustrated on entering into the biaxial phase, apparently due 
to the competing energetic demands of diherent contributions to the Hamiltonian 
arising from incompatible surface anchoring inhuences involving the minor axes. 

The ehect of the variation of the thickness is marginal. The Type - A him seems 
to undergo an anchoring-induced transition in its director structures when the 
is tuned. These two observations could to be useful inputs in designing a practical 
device with switchable in-plane order. 

In summary, we conclude that the Type - A him, involving realistic LC biaxial 
molecular system with simple hybrid anchoring conditions imposed only on the long 
axes is a potential candidate to generate a high degree of in-plane order in the biaxial 
phase, amenable for coupling to external helds. 
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